arXiv:1504.04260v2 [quant-ph] 15 Sep 2015 


Robust quantum correlations in out-of-equilibrium 
matter-light systems 


O. L. Acevedo 

Departamento de Fi'sica, Universidad de los Andes, A.A. 4976, Bogota, Colombia 
E-mail: ol. acevedoSSOuniamdes . edu. co 

L. Quiroga 

Departamento de Fi'sica, Universidad de los Andes, A.A. 4976, Bogota, Colombia 

F. J. Rodriguez 

Departamento de Fi'sica, Universidad de los Andes, A.A. 4976, Bogota, Colombia 

N. F. Johnson 

Department of Physics, University of Miami, Coral Gables, Miami, FL 33124, USA 

Abstract. High precision macroscopic quantum control in interacting light-matter 
systems remains a significant goal toward novel information processing and ultra- 
precise metrology. We show that the out-of-equilibrium behavior of a paradigmatic 
light-matter system (Dicke model) reveals two successive stages of enhanced quantum 
correlations beyond the traditional schemes of near-adiabatic and sudden quenches. 
The first stage features magnification of matter-only and light-only entanglement and 
squeezing due to effective non-linear self-interactions. The second stage results from 
a highly entangled light-matter state, with enhanced superradiance and signatures 
of chaotic and highly quantum states. We show that these new effects scale up 
consistently with matter system size, and are reliable even in dissipative environments. 
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Many-body quantum dynamics are at the core of many natural and technological 
phenomena, from understanding of superconductivity or magnetism, to applications in 
quantum information processing as in adiabatic quantum computing [1], Critical phe¬ 
nomena, defect formation, symmetry breaking, hnite-size scaling are all aspects that 
emerge from the collective properties of the system [2], Spin networks, many-body 
systems composed of the simplest quantum unit, are an obvious starting point to un¬ 
derstand those phenomena, as they enclose much of their complex behavior in a highly 
controllable and tractable way. However, if the system under investigation includes a 
radiation subsystem, new opportunities arise for monitoring and characterizing the re¬ 
sulting collective phenomena [3,4], By devising driving protocols of the light-matter 
interaction, high precision macroscopic control then becomes a possibility, regardless of 
whether the focus is on the matter subsystem, the light, or the composite manipulation 
of both. This is particularly true for the Dicke model (DM) [5], which is the subject of 
the present work. 

The DM describes a radiation-matter system which, despite its simplicity, exhibits 
a wide arrange of complex collective phenomena, many of them specihcally associated 
with the existence of a quantum phase transition (QPT) [6]. Experimental realizations 
of the DM have been presented in different settings, from proposed realizations in cir¬ 
cuit quantum electrodynamics [7], to the recent very successful demonstrations of DM 
superradiance in various cold atom experiments [8]. While the light and matter proper¬ 
ties in the equilibrium ground state are relatively well known [9-13], its fully quantum 
out-of-equilibrium critical behavior is just starting to be understood [14,15]. 

In the DM, both the matter and held are known to act as mediators of an effective 
non-linear self-interaction involving each other [10,13]. These non-linear interactions 
produce interesting phenomena in both atomic and optical systems [16,17]. Among the 
most relevant effects, there is the strong collapse and revival of squeezing [18,19], which 
in many matter states can be related to atom-atom entanglement [20]. Applications 
of such effects are widespread, including high precision quantum metrology [21], and 
quantum information technologies [22]. As the non-linear interactions are only effec¬ 
tive, it is of essential relevance to understand how the eventual matter-held correlations 
could ahect the generation of the desirable quantum squeezing in each subsystem. We 
will theoretically address this problem in the particular setting of the hnite size DM, 
when it is ramped across the critical threshold starting from its initial equilibrium non¬ 
interacting state. 

Some previous works have already examined the dynamical emergence of quantum 
ehects on each subsystem of the DM. However they focused on the semi-classical limit 
and a static coupling after a sudden quench [23-26], small dynamic oscillations around 
a phase space region [14], or under time-delayed feedback control [27]. Our work goes 
beyond this by giving a fully quantum analysis of the process of continuously turning on 
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the interaction, in order to assess the emergent effects of the non-linear self interaction 
in each subsystem. We have carried this out across all dynamical regimes: from the very 
slow adiabatic regime, where the equilibrium ground state results apply, to the sudden 
quench regime which is the initial preparation scheme for most previous non-equilibrium 
results. In between these limits, we have found an unexplored yet very rich intermediate 
set of annealing velocities featuring remarkable amplification of critical quantum effects 
as compared to the near-equilibrium results. This magnification of critical properties is 
followed by a novel chaotic dynamical phase, characterized by giant light-matter entan¬ 
glement [28], negative Wigner quasi-distributions, and signatures of quantum chaos. 

This paper is structured as follows. Section 1 presents the model, dynamical setting, 
and methods of analysis. Section 2 explains the general profile of the dynamical evolution 
under the perspective of matter and light quantum properties, namely squeezing and 
the order parameter (OP). Section 3 gains further insight from the point of view of 
phase space representations of the subsystems. Section 4 establishes the consistency of 
our results, even under dissipative effects, and for different system sizes, by means of a 
power-law relation. Section 5 concludes with some general remarks. 


1. Theoretical framework 


The DM describes a totally symmetric interaction of N matter qubits with a radiation 
mode. Its Hamiltonian is [5], 


H — eJz (jjaf a -(- 2— i^=Jx (a^ T af\ . 
VN 


( 1 ) 


Symbols R = \ Yl!j=i represent collective operators of the qubits, and symbol of 
(a) is the creation (annihilation) operator of the radiation field. Coefficients e and uj 
are single qubit and single mode excitation energies respectively. From now on, we set 
resonant energies, e = a; = 1. The thermodynamic limit (iV —)■ oo) phase boundary 
of the controlled interaction is Ac = \feu:/2 = 0.5, which is slightly different at finite 
values of N [11]. The total system S', represented by a density matrix ps, evolves under 
a unitary part generated by H, and a dissipative part caused by radiation losses of the 
field cavity. The time evolution is then expressed by a master equation [29], 


fL 

dt 


Ps = 


H,ps 


+ K{n + 1) (2dpsd'^ — {afd, /) 5 }) ku (2af psd — {dd\ ps}) , (2) 


where k is the damping rate of the cavity, and n is the thermal mean photon number. 
At initial time f = 0 the system starts at 

m = (3) 

where e“^ = nj (n + 1). Thus, the initial state at A(0) = 0 corresponds to an un¬ 
entangled state formed by the qubit ground state times the field thermal equilibrium 
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state at inverse temperature ft = l/lksT), with ks the Boltzmann constant. The con¬ 
trolled interaction will change with an annealing velocity (AV) r; by a linear ramping: 
X{t) = vt. The interval of interest is A G [0,2], well across the QPT. This work is fo¬ 
cused on subsystem properties. Any subsystem A will be described by a density matrix 
Pa, dehned as the trace with respect to the other degrees of freedom: pA(t) = tr 5 _A {Ps}- 

We explore the dynamical development of the quantum properties in subsystems of 
the matter-light system, caused by this ramped QPT crossing, without resorting to any 
common simplihcations like mean-held, rotating wave or semi-classical approximations. 
Our main results he at AV values outside the well-known adiabatic and sudden quench 
regimes, where quantum correlations such as entanglement and squeezing of each sep¬ 
arate subsystem can get signihcantly enhanced. Exact numerical solutions have been 
obtained by integrating Eq. 2. When k = 0, the evolution lies on the pure Hilbert space 
and is generated by the DM Hamiltonian. Otherwise, the evolution lies in the space of 
density matrices and its generator is the total Liouvillian superoperator of Eq. 2 [29]. 
Thus, the dimension of the numerical evolution vector space is squared as soon as the 
unitary condition is broken. As the dimension of the held Hilbert space is inhnite, the 
Fock basis {|n)} is truncated up to a number where numerical results converge. We 
have taken advantage of every simplifying condition, such as parity conservation when 

K = 0. 


The state of any set of M qubits is the same (all qubits are equivalent) and lies 
on the maximal Dicke manifold = (J -|- 1)J, with J = M/2. The hrst quantum 
property of matter that we address is entanglement, which for a 2 qubit subsystem 
is measured by the Wootters concurrence cw [30]. The second one is spin squeezing, 
which is highly related to multipartite entanglement and usually expressed in terms of a 
squeezing parameter [31]. Under unitary evolution (k = 0), when parity of the total 
system is even and conserved, cw and fq are explicitly related whenever the concurrence 
is diherent from zero [32]: 


l-fl = {N- l)cw 




( 4 ) 


The factor V — 1 in the concurrence is a direct manifestation of the monogamy of 
entanglement [33], since each qubit is equally entangled to any other one and hence 
the hnite amount of possible entanglement is evenly distributed. Any qubit state 
can be visualized by the Agarwal-Wigner function (AWE), which is a Bloch sphere 
representation of the qubit’s density matrix [34], 

N I 

Wq{9,cl>) = J2 

1=0 m=—l 


( 5 ) 
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and Yi^rn are the spherical harmonics. Terms = tr|/)qT;m|, with pg the qnbit 
snbsystem density matrix, are the expected valnes of the mnltipole operator, 

f,„ = ^ (-ijJ-yaTT ( I f, ) \jM) ijM'\ , (6) 

M,M'=-j \ / 

where j = N/2, and ( Tm L m' ) is the Wigner 3j symbol. 

In addition to the matter snbsystems, the DM has the presence of a radiation 
mode that is usnally mnch more experimentally accessible, thanks to the radiation 
leaked by the cavity and properties disclosed by tomographic techniqnes [3]. The 
general state of the mode can be condensed in the form of its Wigner qnasi-probability 
distribntion [35,36], 

OO 

Wb (a, pb) = ^ (-1)'' {n\ (a) pbt) (a) \n ), (7) 

where D (a) = is the displacement operator, and a G C. There is also an 

analogons sqneezing parameter in the held mode, thongh it cannot be directly related 
to a form of entanglement. It is expressed in terms of the variance and covariance of 
held qnadratnres [37], 

= Var{x) + Var{p) — \J {Var{x) — V ar{p))^ + ACov{x,pY. (8) 

The displacement operator in Eq. 7 is related to the qnadratnres by \/2a = x -\- Ip. 

2. Out-of-equilibrium enhancement of quantum correlations 

A general pictnre of the time evolntion of the matter qnantnm properties is revealed by 
means of spin sqneezing (related to matter entanglement) and the expectation valne of 
its OP in Figs. la-b. All the relevant AV are considered in each plot, from the adiabatic 
limit (bottom of vertical axis) to the sndden qnench (top). If the adiabatic condition 
is fnihlled, the spin sqneezing or entanglement of the matter snbsystem corresponds to 
those well known from the eqnilibrinm QPT [9,38]. In the sndden qnench limit, there 
are no appreciable changes since the system essentially stays where it started. In be¬ 
tween those limits, a complex dynamical regime emerges. Near the thermodynamic limit 
eqnilibrinm phase transition A = 0.5, adiabatic evolntion exhibits a maximnm valne of 
entanglement followed by its decay and then the growth of the OP. The anti-correlation 
between sqneezing and OP is present in the dynamical regime as well, so that Figs, 
la and lb are effectively negative images of each other. As the range of high AV is 
entered, the position of the maximnm point of sqneezing and the onset point of the OP 
are pnshed toward higher valnes of A. Any valne of v beyond the adiabatic limit leads 
to the dynamical regime, provided that the controlled interaction is ramped np to a 
high enongh valne. In other words, the sndden qnench condition is only a conseqnence 
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Figure 1. Dynamic profiles (time varies from left to right) of matter subsystem 
properties, (a) Two qubit concurrence cw{N — 1) = 1 — related to spin squeezing 
in Dicke manifold states (see Eq. 4). (b) Scaled order parameter (OP) /A^+ 1/2. 
Results are for N = 81, and unitary evolution (k = 0). Dynamic regimes range from 
adiabatic (bottom) to sudden quench (top) annealing velocities (AV). Compared to 
what is possible in near-equilibrium (lower part), dynamic magnification of squeezing 
occurs at higher AV preceded by the emergence of the order parameter at delayed 
values of A. The insets show the evolution for three selected values of log 2 v (according 
to color) and three different system sizes (according to line style). Strong similarities 
exist for different values of N, pointing toward scaling properties of the results. Roman 
symbols in (b) hint at instants where phase space representations of subsystem states 
are depicted in Fig. 3. 


of the upper bound of A. 

More important than the position of maximum concurrence in the dynamical 
regime, is its value. It reaches up to three times the corresponding value of the adiabatic 
evolution (see Fig. la). This is a remarkable improvement for the collective generation 
of such a distinctive quantum property, and is due entirely to the system being in the 
non-equilibrium dynamical regime. For some given AV values, this magnihed concur¬ 
rence for different number of qubits N follows closely the scaling cw{N — 1) (see inset 
of Fig. la), keeping the spin squeezing parameter curves almost size independent, and 
gets modestly better as the system size grows. Some size dependent deviations are 
present in curves for the log 2 U = —1.56 case. The sudden death point of concurrence 
is also virtually size-independent. The scaling properties of the dynamical regime will 
be discussed in section 4. The matter OP, for some chosen AVs, also shows very good 
scaling collapse in its dynamical evolution, i.e. the curves in the inset of Fig. lb are 
essentially size independent thanks to the scaling Jz/N in the OP. 

With respect to the maximum matter OP in Fig. lb, it is bounded. Therefore 
independent of the regime, it cannot be increased. The kind of enhancement revealed 
by the OP is a magnihcation of superradiance, which is the traditional way to describe 
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the QPT of DM [39]. The OP indicates the scaled number of excitations within the 
subsystem, which is bounded because of hniteness of N. In the dynamical regime, 
excitations develop more quickly (see blue curve in inset of Fig. lb). In the case of 
eventual leakage of excitations through the cavity (see Section 4 for cavity loss effects), 
it manifests itself as a sharp burst of photons which are suddenly released. Bursts like 
these are excellent indicators of the DM QPT in cold atom experimental realizations [8] . 

Qubit quantum properties are very important for applications in quantum compu¬ 
tation and quantum information. However, the hnite character of the matter subsystem 
conceals much of the complexity of the evolution. The held subsystem, with its own 
relevance in highly controlled quantum optics, does not have this restriction. Field 
dynamical prohles analogous to that of Fig. 1 are shown in Figs. 2a-b. Similarities 
with Figs, la-b are noticeable, though with relevant differences. A large dynamical 
magnihcation of the held quadrature squeezing is also present, but it dies oh later than 
its matter counterpart. Thus the development of the held OP is delayed as compared 
to that of the qubits. Superradiance, as seen by value of the held OP, is now not only 
sharper, but it can be as much as twice the value attainable with adiabatic ramping at a 
given A. This greatly enhances the intensity of superradiant bursts. Radiation intensity 
enhancements of this kind have also been predicted in the DM submitted to sudden 
quenches and a lossy cavity [36]. Dynamically enhanced and suppressed superradiance 
alternate after the hrst burst (see blue curve in inset of Fig. 2b), a behavior absent in 
the equilibrium QPT. A wavy plateau in hgure lb coincides with that alternating stage. 
The dynamical phase marked by these oscillations is characterized by giant light-matter 
entanglement, with emergent quantum properties beyond the critical ones [28]. 

Magnihcation of squeezing and the other ehects discussed so far are caused by a 
prolonged realization of ehective non-linear interactions within both subsystems, after a 
preparation stage where the sudden quench approximation holds. Then, the non-linear 
processes are realized as one axis squeezing with a transverse held for matter [19,31], 
and Kerr-like interactions for radiation [17]. During these separate squeezing processes, 
both radiation and matter act as ehective interaction mediators for each other. This 
extends to the non-equilibrium case static correspondences with matter-only [10], or 
radiation-only systems [13]. Separate squeezing, which is related to internal quantum 
entanglement, suddenly dies as soon as combined matter-light entanglement emerges, 
because of monogamy and the breaking of the ehective non-linear interaction condition. 
As the held has more information capacity, its squeezing can survive longer. By con¬ 
trast, in the adiabatic regime, entanglement within and between subsystems compete 
against each other, because they occur simultaneously at the critical point. 
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Figure 2. Corresponding dynamic profiles and insets of Figs, la and 2b, now for the 
field subsystem, (a) Evolution of 1 — as defined in Eq. 8, whenever it is greater 
than zero (squeezed radiation), (b) Scaled OP of light (a’^a) /N. Major tendencies 
of Fig. 1 are replicated here. However the unbounded nature of the field OP allows 
an appreciation of the field intensity’s oscillation (associated with supperradiance) 
around the equilibrium values, with intervals of significant dynamical enhancement of 
superradiant bursts. The oscillatory behavior is related to a combined light-matter 
quantum chaotic stage. 


3. Matter-light quasi-probability behavior 

Complete state representations of the subsystems state reveal more details of the pro¬ 
cesses involved in the dynamical regime. Figures 3 and 4 show several snapshots of 
Wigner and AWF quasi-distributions at different instants. The enormous difference in 
the squeezing amount between the dynamical regime and the adiabatic one, as well as 
a difference in the direction of squeezing, is now graphically clear. This can be seen 
by comparing Figs. 3-4 (iv) with Figs. 3-4 (ii), and was previously shown for differ¬ 
ent AV values in Figs.l and 2. In any regime, death of squeezing is associated with 
a splitting in half of the distributions, caused by the spontaneous symmetry breaking 
of the QPT [28]. This splitting is symmetrical along the x and —x direction in both 
subsystems. The system’s density of excitations, indicated by the OP, increases as the 
distributions become displaced away from the initial state. 

Features of the complex light-matter entangled stage are clarihed by Figs. 3-4 (v- 
vi), which contrast with the adiabatic splitting in Figs. 3-4 (hi). The dynamical splitting 
of the distribution leads to negative scars and complex patterns for both subsystems, 
though the phenomenon is far more conspicuous in radiation. Donut shapes with a neg¬ 
ative centered AWF such as that in Fig. 3 (vi), have been experimentally obtained in 
3000 atoms with just a single photon [4], In addition, round-tailed interference patterns 
such as that in Fig. 4 (vi) have been obtained for light in Kerr-like media following 
a Fokker-Planck equation [18], which confirms the presence of non-linear effective in- 
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Figure 3. Agarwal-Wigner Wq at the instants indicated in Fig. lb by the 
corresponding number (dissipative effects are ignored, k = 0). They are phase space 
representations of the matter subsystem, depicted on the Bloch sphere. Color scale of 
all Wigner functions can be very different. Maximal dynamical squeezing in (iv) is 
much greater than the adiabatic counterpart in (ii). Also, negative complex structures 
develop in the combined light-matter chaotic stage in (v-vi); which contrast with an 
adiabatic ordered state as in (iii). The stage in (iii,v,vi) has also no squeezing and 
no concurrence, since the Wigner function is now split along the x and —x directions 
and no longer concentrated around the initial state. Opposite Bloch hemispheres are 
not shown because of symmetry: Wq{9, cj)n) = Wq{9, cj)). 


teractions. These results mean that the field is able to exhibit chaotic behavior itself, 
regardless of its entanglement to the matter subsystem. Despite being a single radiation 
mode, the field acts as a reservoir that dissipates the quantum correlations present in the 
squeezed states of each subsystem. If the qubits were not coupled to the field, squeezing 
could have revivals after its sudden death [19]. The field Wigner function is not only 
full of negatives regions (a marker of non-classical light), but it also contains abundant 
so-called sub-Planckian structures which have been related to quantum chaos [40] . The 
finite-size DM is non-integrable, and its ordered phase has been connected to chaotic 
features [6,41]. This chaotic onset is responsible for the small size dependence on curves 
in insets of Figs. 1 and 2 for the log 2 U = —1.56 case, since the symmetry breaking 
occurs deeper in the ordered phase. As spin squeezing is very sensitive to chaos [26], it 
is no surprise that its sudden death becomes irreversible once the light-matter entangled 
stage arises. 
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Figure 4. Analogous pictures of Fig. 3 but for the field Wigner distribution Wi,. This 
is represented in the x — p plane of position (vertical) and momentum (horizontal) 
quadrature. All figures have the same x — p scale but color scale can be very different. 
The chaotic stage in (iv-v) exhibits sub-Planck structures. Negative regions of x are 
not shown because of symmetry: Wb{x,p) = Wb{—x,—p). 


4. Scaling and robnstness of non-equilibrinm quantum correlations 
4.1. Robustness of results 

So far, unitary evolution have been assumed in this work. Figures 5a-e address how 
the presence of losses in the cavity affect the dynamics. All the main results found in 
K = 0 evolution survive very well if decoherence towards the environment is two or¬ 
der of magnitudes lower than the main energy scale. Even if dissipation is at values 
of just an order of magnitude below, spin squeezing effects remain highly robust, with 
increasing noise resistance with system size. Field squeezing surprisingly survives to 
dissipation regimes comparable to the Hamiltonian dynamics itself. On the other hand, 
detailed features of the chaotic stage (such as OP oscillations, negative regions, and 
sub-Planckian structures) are far more sensitive to decoherence, requiring losses to be 
at levels of k = 0.01. These very sensitive features have been proposed as tools for mea¬ 
suring very weak forces [42] . In our analysis, we have found that small finite values of n 
(such as those typical at the ultra-low temperatures of most experimental realizations) 
do not change dramatically the conclusions; they just slightly intensify the process of 
decoherence. Previous works on non-equilibrium DM under different settings have also 
confirmed that some major phenomena present here, such as spin squeezing, survive 
losses. Enhanced superradiance can be detected [36], and squeezing can still have time 
to develop [18,24], We note that, despite perturbing the quantum state of the system, 
a lossy cavity compensates by providing the experimental possibility to monitor the 
dynamic evolution of the field in time [3]. 
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Figure 5. Effect of a leaky field cavity (expressed in terms of cavity decay rate k) 
on the dynamical enhancement of quantum properties, (a) Wigner distribution for 
AV log 2 V = —1.58, = 15 qubits at instant A = 2. (b-e) Different field and matter 

properties as function of time for log 2 v = —3.86, and three different system sizes. 
Quantum chaos signatures (sub-Planckian structures, negative Wigner distributions, 
and oscillations in OPs) are highly sensitive to dissipative perturbations. On the other 
hand, squeezing and internal entanglement effects are far more robust, especially for 
field squeezing. In all figures fi = 0. 


The robustness of quantum correlations suggests that some of the effects described 
in this paper may be accessible under current experimental realizations of the DM. In 
particular, Klinder et ah have proven that the kind of control described here as X{t) 
can be realized as a pump laser power increasing over time, and that the noise present 
in current experiments (nearly analogous to our k = 0.1 case) is small enough that the 
squeezing related quantum correlations predicted in our work could be experimentally 
assessed [43]. As held squeezing is the most robust quantum correlation (see Fig. 5e) 
and it is experimentally accessible thanks to a leaky mirror, tomographic techniques 
and photon statistics could be the hrst front for experimental comparisons with our 
predictions [3]. Other tomographic techniques could also analyze properties of the 
matter subsystem state [44], Given the specihc degrees of freedom forming the qubit 
states in current realizations (basically matter wave modes) [43], suitable equivalences 
of those state resolving techniques to this kind of degrees of freedom should be devised. 

4-2. Scaling of results 

Our results exhibit substantial size invariance, as the insets of Figs. 1 and 2 show. Even 
the small system sizes in lossy cavity results in Figs. 5b-e are highly size independent. 
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which justifies our claim that the conclusions drawn there can be taken as general 
trends. Figures 6a-b show that the onset point of the dynamical regime occurs at points 
\d following a power law, (A^ — Ac) oc This kind of scaling relation is typical in 

QPT critical properties [2], This scaling laws extends our previous results on adiabatic 
scaling theory [15] to the fast regime. Even though we have used the same kind of scaling 
variables as in the near adiabatic case, many results in the dynamic regime collapse very 
well without any re-scaling of A. 



Figure 6. Size-independence of the instant Ad when field and qubits organize in a 
combined light-matter excited state, (a) Instant of a sudden death of squeezing, in both 
matter and radiation subsystems, (b) Instant when each system reaches a significant 
value of its OP. These values are /-^ + 1/2 = 0.1 and {a)a) /N = 0.0123. At 
sufficiently high AVs, all the instances chosen are related by the same scaling power 
law relation (Ad — Ac) oc 

The origin of the scaling exponent can be understood in a rather simple way, using 
concepts from the Kibble-Zurek mechanism [1]. The scaled time is the instant where 
the healing time of the system (as measured by the inverse mean energy gap in the 
spectrum 5“^) has the same order of magnitude, i.e. when M ~ 1. As 5 oc (A^ — Ac)^^^ in 
the ordered phase [10], and t oc (A^ —Ac)n“^, we get the exponent 2/3 present in all lines 
of Figs. 6a-b. The fact that equilibrium critical exponents have been used reveals that, 
even though the emergence of the dynamical regime is well inside the A > 0.5 range, its 
mechanism is still related to the QPT, and so many critical phenomena get dynamically 
enhanced. On the other hand, the fact that the evolutions are very similar in the chaotic 
region, despite the different sizes, has to do with the relatively equal structure of the 
interaction well inside the superradiant phase, provided it has not bounced against the 
finite limits of the matter subsystem. 
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We investigated the details of the microscopic properties of both matter and light sub¬ 
systems of the DM when macroscopically driven into the superradiant phase from an 
initial non-interacting equilibrium state. We have found that, in between the traditional 
driving schemes (either near-adiabatic or sudden quench), a previously unnoticed inter¬ 
mediate annealing regime emerges, and that this regime can be seen as advantageous 
for many quantum control processes. The onset of this intermediate phase is marked 
by induced non-linear self-interactions in both matter and held subsystems, developing 
squeezed states in each of them that are related to internal entanglement. This squeez¬ 
ing process is much stronger than what can be achieved in near-adiabatic evolutions, 
since the internal entanglement does not have to compete against the matter-held en¬ 
tanglement that arises in later stages. 

We have also found that the dynamical squeezing mechanism is succeeded by 
the development of a chaotic and entangled radiation-matter dynamical phase. Clear 
signatures of quantum chaos such as fragmented Wigner functions, have been identihed 
in this phase. This stage has its own dynamically enhanced properties such as 
superradiance, phase order, and sensibility to weak forces. We have established the 
robustness of all the quantum enhancing processes under dispersive conditions and its 
invariance to system size. We have shown that this last invariance condenses into a 
power-law relation associated with critical exponents. We hope this understanding will 
prove important for developing schemes to generate squeezed and other entangled states, 
as have been proposed in the contexts of quantum metrology and quantum information 
processing. 
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